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ABSTRACT 

Neutrino-cooled hyperaccretion disks around stellar mass black holes are plau- 
sible candidates for the central engine of gamma-ray bursts. We calculate the 
one-dimensional structure and the annihilation luminosity of such disks. The neu- 
trino optical depth is of crucial importance in determining the neutrino cooling 
rate and is in turn dependent on the electron fraction, the free nucleon fraction, 
and the electron degeneracy, with given density and temperature of the disk 
matter. We construct a bridging formula for the electron fraction that works 
for various neutrino optical depths, and give exact definitions for the free proton 
fraction and free neutron fraction. We show that the electron degeneracy has 
important effects in the sense that it enlarges the absorption optical depth for 
neutrinos, and it along with the neutronization processes favored by high tem- 
perature cause the electron fraction to drop to be below 0.1 in the inner region of 
the disk. The resulting neutrino annihilation luminosity is considerably reduced 
comparing with that obtained in previous works where the electron degeneracy 
was not considered and the electron fraction was simply taken to be 0.5, but it 
is still likely to be adequate for gamma-ray bursts, and it is ejected mainly from 
the inner region of the disk and has an anisotropic distribution. 

Subject headings: accretion, accretion disks - black hole physics - gamma rays: 
bursts - neutrinos 



1. Introduction 

Theoretical models for gamma-ray bursts (GRBs) can be divided into two categories: 
those named fireball models that treat the shock interaction of relativistic outflows and 
production of gamma rays and afterglows in other wavelengths (see, e.g., Meszaros 2002 and 
Zhang & Meszaros 2004 for reviews), and those that explore the central engine of the fireball. 
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i.e., the energy source of relativistic outflows. Most popular models in the latter category 
are in common invoking a stellar-mass black hole accreting with a hypercritical rate, of the 
order of 1 Mq s~^. The main problem in these models is how to convert some fraction of the 
released gravitational energy of the accreted matter into a relativistic outflow, creating an 
explosion with energy ~ 10^''-10^^ ergs (depending on whether emission is isotropic or it is 
beamed). Two mechanisms have been proposed to tackle this problem: the neutrino emission 
and annihilation, and the energy extraction from the accretion disk and/or the black hole 
via magnetohydrodynamical processes (see, e.g., Popham et al. 1999 and Di Matteo et al. 
2002 for references). The former mechanism is easier to understand and can be calculated 
more accurately, and is the topic we wish to discuss in this paper. 

In the inner region of such a hyperaccretion disk the density and temperature are so 
high (p ~ 10^° g cm~^, T ~ 10^° K) that photons are totally trapped, and large amounts of 
energetic neutrinos are emitted, carrying away the viscous dissipation energy of accreted gas. 
Annihilation of some fraction of emitted neutrinos produces a relativistic electron-positron 
outflow to power a GRB. The properties of such a neutrino-cooled accretion disk model 
were first worked out in details by Popham et al. (1999). From the observational point of 
view, the key question to be answered in this model is whether the neutrino annihilation can 
indeed provide sufficient energy for GRBs. Popham et al. (1999) gave a positive answer to 
this question, but they assumed a priori that the accretion disk is transparent for neutrinos, 
thus their neutrino radiation luminosity (before annihilation) Lj, and accordingly neutrino 
annihilation luminosity L^^ might be overestimated. Later, Di Matteo et al. (2002) recal- 
culated Li, and Li^x? in the neutrino-cooled accretion disk model, taking the neutrino opacity 
into account. They obtained that even for a modest mass accretion rate the inner region 
of the disk becomes opaque and neutrinos are sufficiently trapped, the resulting L,^j7 is < 
10^° ergs s^^ and is inadequate for GRBs. In a recent work (Gu et al. 2006) we showed 
that when the general relativistic effect is considered and the contribution from the optically 
thick region is included, the neutrino-cooled accretion disk can work as the central engine 
for GRBs, although the correct L,yTy is somewhat lower than that of Popham et al. (1999) 
because of the effect of neutrino opacity. 

In addition to the general relativity and the neutrino opacity, there are certainly other 
factors that may influence the neutrino radiation and annihilation of a neutrino-cooled ac- 
cretion disk, such as the electron degeneracy and the electron fraction. In our previous work 
(Gu et al. 2006) the treatment of electron degeneracy was oversimphfled, and the electron 
fraction was simply taken to be 0.5, i.e., assuming an equal mix of protons and neutrons 
since ~ Up/ (rip -|- rin) , where rip and nn are the total number density of protons and that 
of neutrons, respectively. Kohri & Mineshige (2002) pointed out that when electrons are 
degenerate, there emerges an important consequence that the electron-positron pair creation 
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and accordingly the neutrino emission are suppressed. Kohri et al. (2005) took great care 
to calculate the ratio of free neutrons to free protons h-^Jn^ since this ratio has a large effect 
on the neutrino emission rates: the true ratio n^lh^> 1 (or < 0.5) resulted from neutron- 
ization processes will also lead to a suppression of neutrino emission (note that they did not 
distinguish between Un, Up and §2.3.1). They made calculations for the neutrino 

emission rates and other quantities even in the delicate regime where the electron degeneracy 
is moderate, which is also a significant improvement over previous works. Lee et al. (2005) 
also considered the effects of electron degeneracy and electron fraction in neutrino-cooled 
accretion disks: they used an expression for the pressure of ultra-relativistic electrons with 
arbitrary degeneracy and calculated with an approximate bridging formula that allows for 
a transition from the neutrino optically thin to optically thick regime. These works (Kohri 
& Mineshige 2002; Kohri et al. 2005; Lee et al. 2005) made advances in microphysics, but 
they were within the Newtonian framework. Very recently, Chen & Beloborodov (2006) 
presented fully relativistic calculations of the structure of neutrino-cooled accretion disks 
around Kerr black holes and proved that both the electron degeneracy and the electron 
fraction dramatically affect the disk and its neutrino emission. 

It is seen from the above brief review that the electron degeneracy and the lower elec- 
tron fraction certainly suppress the neutrino emission and reduce the neutrino annihilation 
luminosity of a neutrino-cooled accretion disk. The question remains whether the reduced 
neutrino annihilation luminosity is still adequate for GRBs. In this paper we try to refine 
our previous results of the structure and luminosity of neutrino-cooled accretion disks by 
considering the relevant microphysics more completely and more accurately. 

2. Physics of Neutrino-cooled Accretion Disks 

2.1. Hydrodynamics 

We limit the central accreting black hole to be a non-rotating one, its general relativistic 
effect is simulated by the well-known Paczyhski & Wiita (1980) potential $ = —GM/{R — 
Rg), where M is the black hole mass, R is the radius, and Rg = 2GM/c^ is the Schwarzschild 
radius. 

As in the relevant previous works, the hydrodynamics of hypcraccrction disks is expected 
to be similar to that of normal accretion disks in X-ray binaries, and is well approximated 
by that of steady axisymmetric height-averaged accretion flows (e.g.. Chap. 3 of Kato et al. 
1998). Accretion in the disk is driven by viscous stress, and the kinematic viscosity coefficient 
is expressed as v = ac^H, where H is the half thickness of the disk; Cg = {P/ pY^"^ is the 
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isothermal sound speed, with P and p being the pressure and mass density, respectively; 
and a is a dimensionless constant parameter that absorbs all the detailed microphysics 
of viscous processes. The angular velocity is approximately Keplerian, i.e., VL = f^K = 
{GM / R)^/"^ / [R — Rg). The disk is in the vertical hydrostatic equilibrium, and this gives 
H = Cs/^Ik- With these simplifications the problem is reduced to be one-dimensional, i.e., 
all physical quantities depend on R only. 

The constant mass accretion rate M is expressed from the continuity equation as 

M = -4:7rRHpv, (1) 
where v is the radial velocity that can be read from the angular momentum equation as 

v = aHc,{l-j/n^RY'^^, (2) 

where j is an integration constant determined by the zero-torque boundary condition at the 
last stable orbit, and it represents the specific angular momentum (per unit mass) of the 
matter accreted into the black hole. 



2.2. Thermodynamics 

The energy equation is generally written as the balance between the viscous heating 
and the cooling rates (per unit area of a half-disk above or below the equator), 

Qvis = Q-. (3) 

The viscous heating rate Qvis is similar to that of normal accretion disks, 

Qvi, = -^Mn^R{i-j/n^R^)^. (4) 

47r an 
But the cooling rate Q~ is crucially different, it has three contributions: 

Q — Qphotodis + Qadv + Qu- (5) 

In this equation there is no cooling term of photon radiation (it is practically zero in our 
calculations). Instead, photons are totally trapped in the disk, so they contribute to the 
advective cooling and the pressure (see below). The coohng rate by photodisintegration of 
a-particles Qphotodis is 

Qphotodis = 6.8 X lO^^piovH ^^^" cgs units, (6) 
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where pio = p/lO^^g cm ^, and Xnuc is the mass fraction of free nucleons (e.g., Kohri et al. 
2005). The advective coohng rate Qadv is 

where s is the specific entropy, ^ oc — dlns/dlni? is taken to be equal to 1, i.e., ds/dR is 
approximated as s/R (Kohri & Mineshige 2002), and is the mean mass of a nucleon. 
The entropy of degenerate particles is small and can be neglected. The three terms in the 
brackets of equation (7) are the entropy density of photons, of free nucleons and a-particles, 
and of neutrinos, respectively; and is the energy density of neutrinos, for which we adopt 
a bridging formula valid in both the optically thin and thick regimes (Popham & Narayan 
1995; Di Matteo et al. 2002), 

^ _ ^ (7/8)aT\rj2 + 1/^3) 

where Tj,^ is the total optical depth for neutrinos, Ta^^,. is the absorption optical depth for 
neutrinos, and the subscript i runs for the three species of neutrinos , , and Ut- . The 
cooling rate due to neutrino loss Qi, is expressed in accordance with the above equation (Kohri 
et al. 2005), 

Q =V (9) 

V(3/4)[r../2 + l/y3 + l/(3v.J]- 

The equation of state is also very different from that of normal accretion disks, as 
the contributions to the pressure from degenerate electrons (we assume that nucleons are 
not degenerate throughout the present paper) and from neutrinos should be included, it is 
written as 

P = Pgas + i'rad + Pc + P.- (10) 

The gas pressure from free nucleons and a-particles Pgas is 

_ kspT 1 + 3Xnuc , . 

m„ 4 ■ ^ ' 

The photon radiation pressure Prad is 

Pr.d = aPyS. (12) 



The electron pressure is from both electrons and positrons and should be calculated 
using the exact Fermi-Dirac distribution. No asymptotic expansions are valid because at 
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different radii electrons may be with different degrees of degeneracy and may be relativistic 
or nonrelativistic. It reads 

Pe = A-+A+, (13) 

with 



1 /-oo 4 I 
Pe^ = TT^TUT^ / dp , ^ , , (14) 

where fic is the chemical potential of electrons, and the electron degeneracy is measured by 
the degeneracy parameter defined as ri^ = j-ic/k-QT . We agree with Lee et al. (2005) that since 
the presence of relativistic e^e^ pairs is automatically taken into account in the expression 
for Pe, there is no alteration to the numerical factor 1/3 in the expression for Prad, nor to 
the factor 4/3 for the photon entropy in equation (7). The neutrino pressure is 

Pu = uj2,. (15) 



2.3. Microphysics 

All the equations in the above two subsections can be combined into only two equations, 
i.e., equations (1) and (3). In these two equations there are seven unknown quantities, namely 
p, T, X-aaci 'Tvii^i — ^c^^ht^t), and /ie(or ?7e). Therefore, one has to find more equations 
relating these unknowns with the knowledge of microphysics. 

2.3.1. Neutrino optical depth 
The total optical depth for neutrinos is 

Ti^i = rs,i., + Ta,i.,. (16) 

The optical depth for neutrinos through scattering off free nucleons, a-particles, and 
electrons Tg^ut is given by 

Ts,ui = H/Xi,^ = H[ap^^.np + an,uinn + cra,vina + cre,n{n^- + rie+)], (17) 

where X^^ is the mean free path; ap^y^ , (Jn,ui , o'a,Ui and (Xe^ui are the cross sections of scattering 
on protons, neutrons, a-particles, and electrons; rip , rij^ , Ua , n^- and 1%^+ are the number 
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densities of free protons, free neutrons, a-particles, electrons, and positrons, respectively. 
The four cross sections are given by (Burrows & Thompson 2002) 



4 



[{Cv,..-l? + ^9\{CA,.,-lfl (18) 



^ (19) 



a^^^^^Aa^ElsiuHw, (20) 



where ctq = 1.76 x 10~^^cm^,£'^,. is the mean energy of neutrinos in units of (mgC^), qa ~ 1-26, 
sin^ Ow ~ 0.23, Cv,.^ = 1/2 + 2 sin^ Qw , Cy,,^ = Cy,.^ = -1/2 + 2sm2% , C^,,^ = C^.p, = 
CA,v-r = 1/2 , CA,Fe = CA,Uf, = CA,ur = -1/2- Slnce Up = np - 2na and = Un - 2na, 
the free proton fraction 1^, the free neutron fraction Y^, and the a-particle fraction are 



related to Y^ and X^uc ^ 



Yp = = n - (22) 



i; = ^^ = l-n-^-^, (23) 
rip + rin 2 



= = (24) 

rip + nn 4 



and ne- , ne+ are given by the Fermi-Dirac integration, 

I , o 1 



n^T = I dpp^ , (25) 

(cf. Kohri et al. 2005, they took -hp = Up and hn = n^, and thus Yp — Y^ and Y^ — 1 — Y^, 
which are valid only for the completely dissociated matter, i.e., X^uc — 1)- 
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The absorption depth for neutrinos Ta,ui is defined by 

^-'•^^ = 4(7/8)aT4' ^^^^ 
where g^,. is the total neutrino coohng rate (per unit volume) and is the sum of four terms, 

(27) 

The neutrino cooling rate due to the URCA processes guRCA relates only to and is repre- 
sented by the sum of three terms (Chap. 11 of Shapiro & Teukolsky 1983; Yuan 2005), 

(28) 

with 

q2 QQg2 poo 

qp+e-^n+.. = ^2^(1 + ^9l)np dEe Ee ^/E,^-mM{Ee - Qffe-, (29) 

Q'2 poo2 n poo 

g„+e+^p+Pe = (1 + ^dDfin J ^ dE, E, VEe'-m^^c^iEe + Qffe+ , (30) 

qn^,+e-+v^ = ^^^(1 + ^gDnnj"^ ^ dE, E, ^/E,^-mM{Q - E,)\l - /e-), (31) 

where Gp ~ 1-436 x 10~^^ergs cm^, cos^ 9c ^ 0.947, Q = {m^-m^)(?, and /e=F = {ex?^[{E^^ 
He)/kBT] + l}"^is the Fermi-Dirac function. The last term in the right hand side of equation 
(28), i.e., equation (31), is small comparing with the other two terms, and was usually not 
included in the literature. The electron-positron pair annihilation rate into neutrinos ^6-6+ 
is (e.g., Itoh et al. 1989) 

ge-e+-..eFe ~ 3.4 X 10=^3^1^1 ergs cm-=^ s-\ (32) 

qe-e+^.,v, = qe-e+^.^vr ~ 0.7 X lO^^^y^^ crgs cm-3 s-\ (33) 

where Tn = T/IO^^K. Expressions (32) and (33) are valid in the nondegenerate limit 
r/c <^ 1, and when electrons are degenerate gc-c+ becomes negligible. The nucleon-nucleon 
bremsstrahlung rate gbrcm through the processes n + n ^ n + n + Ui + Vi is the same for the 
three species of neutrinos (Hannestad & Raffelt 1998; Burrows et al. 2000), 

gbrem ~ 1-5 X lO'V^o^u" ergs cm-=^ 8"^ (34) 
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As to the plasmon decay rate Q'piasmon, only that through the process 7 — > + needs to 
be considered, where plasmons 7 are photons interacting with electrons, 

7^ 

gpia^mon ~ 1-5 X lO^^T^^^^ exp{-^p){l+^p){2 + -^) crgs cm-3 s-\ (35) 

where 7p = 5.565 x 10~^[(7r^ + 3r7g)/3]^''^(Ruffert et al. 1996). It is expected that gfbrem and 
?piasmon Can bccome important only at very high electron degeneracy. 



2.3.2. Electron Fraction 

It is seen that the electron fraction 1^ mentioned in Introduction has appeared in the 
expression of neutrino optical depth T„^ . In order to calculate r^,. and accordingly the neutrino 
coohng rate (see eq. [9]), more knowledge about Ye is required. 

Beloborodov (2003) proved that /3-cquilibrium among free neutrons, free protons, and 
electrons is established in disks with M > 1O^^(q;/O.1)^^^(M/M0)^/^ g s^^, a condition that 
is hkely to be fulfilled for hyperaccretion disks; and noted that a distinction needs to be 
made to determine the equilibrium composition depending on the optical depth of the disk 
material. Kohri et al.(2005) discussed various timescales in neutrino-cooled accretion disks 
in details, and showed that for sufficiently large M and not too large R the timescale for the 
reactions from proton to neutron and from neutron to proton is shorter than the dynamical 
timescale (the accretion time), i.e., the /^-equilibrium is likely to be realized. 

If the disk material is opaque to neutrinos, there are reversible reactions e~ +p ^ n + i/e 
and e^ + n ^ p + V^. The chemical potential of neutrinos can be ignored because the number 
density of neutrinos and that of antineutrinos are likely to be equal, then the /^-equilibrium 
condition is 

Afn = A*p + fJ'e, (36) 

where /i^ and /ip are chemical potentials of neutrons and protons, respectively. On the other 
hand, if the material is transparent to neutrinos, there are no reversible reactions as for the 
neutrino opaque material, but the /^-equilibrium can be reached when the rate of reaction 
e~ -|- p — > n -|- z/c is equal to that of reaction e'^ + n ^ p + 1/^. Yuan(2005) calculated these 
two rates and obtained 

A«n = A^p + 2//e- (37) 
The neutron-to-proton ratio in /3-equilibrium is given by 

■z^ = exp [(/^n - A*p - Q)/kBT] , (38) 

Up 
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which results in 



"p 

for the neutrino opaque limit, and 

lg?. = %^, (40) 

for the neutrino transparent limit. In order to allow for a transition from the optically thin 
to optically thick regime, we adopt a treatment similar to that in Lee et al. (2005), i.e. 
introducing a weight factor /(rjy) = exp(— r,yg) and writing in a combined form, 

Ig?^ = + [1 - /(r.)]^^. (41) 

By using the relation — n-p/{n-p + nn) and equations (22), (23), and (24) we finally arrive 
at 

Fe = -(1 - ^nuc) + i^^^^^r^^ 7^- (42) 

Some comments on equation (42) are in order. First, this equation is more exact than 
equation (12) of Lee et al. (2005), theirs is the first term of the expansion of ours. Second, 
we notice that Kohri et al. (2005) discussed the ratio fi^jfip at length. They started with 
considering the transition rates from proton to neutron and from neutron to proton, then 
obtained the expression of n-ajn^ , i.e., equation (39) for the neutrino opaque limit, and 
adopted an approximate procedure to estimate this ratio for the neutrino optically thin 
regime. We stress that equation (40) is rigorous for the neutrino transparent limit, since it 
is also derived from the transition rates between protons and neutrons (Yuan 2005); then 
we construct equation (42) as a bridging formula between the opaque and transparent limits 
and expect it to work for the regime where the neutrino optical depth is moderate. 

An additional relation between and Xnuc can be found from the equation of nuclear 
statistical equilibrium (e.g., Meyer 1994), 

Y{Z,A) = G{Z,A%{2;)^-\'^^-'''^IH''''-''^l^\A''l\^ 

where Y(Z,A) is the mass fraction of a kind of particles with the charge number Z and 
mass number A, G{Z, A) is the nuclear partition function, C(3) is the Riemann zeta func- 
tion of argument 3, = [2({3){kBT)^]/[K'^{kc)^pNA\ is the photon-to-baryon ratio, Na is 
the Avagadro constant, and B{Z,A) is the binding energy of the nucleus. As in all the 
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previous works, we assume that all heavy nuclei are a-particles. This should be a reasonable 
assumption because all nuclei heavier than a-particles contain approximately equal numbers 
of neutrons and protons. Then by using equations (22), (23), and (24), equation (43) yields 

4(n - ^^^)(l - n - ^^^)(l - ^nuc)-'/' = 8.71 X lOVuf^^T^f exp (-1.64/Tn).(44) 

Only Chen & Beloborodov (2006) adopted this equation, while other previous works used 
a simple expression of X^uc that was obtained by taking 1^ = 0.5 and did not reflect the 
interdependence between 1^ and X^^c (e.g., eq. [40] of Kohri et al. 2005). 



2.3.3. Electron Chemical Potential 

The electron chemical potential is determined by the condition of charge neutrality 
among protons, electrons, and positrons, 

np = — ^ = ne--ne+, (45) 
rriu 

with Ue- and ne+ given in equation (25). 



2.4. Summary of Equations 

The system of equations is closed, as there are eight equations, i.e., equations (1), (3), 
(42), (44), (45), and three equations (16) for eight independent unknowns p,T,ri,^,rt,^,r,^^,Xnuc, 
Ye, and //e(or rje), which can be numerically solved as functions of R with given constant 
parameters M.M.a, and j. All the other quantities such as P and its components, Q~ and 
its components, and composition fractions Yp, 1^, and are obtained accordingly. 

Our equation set has the following advantages: (1) Effects of relevant factors are taken 
into account in a combined way, including the general relativity, the inner boundary condition 
of the disk, various processes that contribute to the neutrino cooling and the neutrino opacity, 
the electron degeneracy, the electron fraction, and the coexistence of electrons, positrons, free 
protons, free neutrons, and a-particles. (2) Whenever the electron degeneracy is concerned 
we use the exact Fermi-Dirac integration rather than the analytical approximations that 
are valid only for extreme cases. (3) We take great care to calculate the neutrino optical 
depth and the electron fraction. In doing so, we make a careful distinction between the total 
nucleon number densities n^, Up and the free nucleon number densities nn,np, so that the 
composition of the disk matter can be exactly described by fractions , Yp , , and Y^ ; 
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we propose a new bridging formula for Ye (eq. [42]) from the /3-equilibrium condition, which 
is apphcable to both the neutrino optically thin and optically thick regimes; and we adopt 
equation (44) that connects with X^uc- 



3. Numerical Results for The Disk Structure 

This section presents our results for the structure of a neutrino-cooled accretion disk, 
obtained by numerically solving the set of eight equations in §2. Figures 1-6 show physical 
quantities of the disk matter as functions of R. In all these figures the necessary constant 
parameters are fixed to be their most typical values, i.e., M = 3Mq, M= 1 Mq s^^ (except 
for Fig. 5 where M= 5Mq s~^),a = 0.1, and j = l.ScRg (see, e.g., Gu et al. 2006 for the 
discussion of j) . 

Figure 1 is for the density p, temperature T, and electron degeneracy parameter 1]^. 
It is seen that from R = 500Rg inward to R = 3Rg, p increases by about four orders of 
magnitudes (Fig. la), T increases by about eight times (Fig. lb); and in the innermost 
region of the disk p reaches to ~ 10^^ g cm"'^, and T reaches to ~ 3 — 4 x 10^°K. With such 
densities and temperatures, r/e increases first with decreasing i?, reaches to its maximum 
value at i? ~ 65it!g, and then decreases with decreasing R because of increasing temperature 
(Fig. Ic). The behavior of rfe obtained by us is consistent with that of Kohri et al. (2005, 
Fig. 4 there); but is not so with that of Chen & Beloborodov (2006, Fig. 3 there), who 
got that 77e always increases with decreasing R in the very inner region. We think that r^e 
ought to decrease when temperature is very high. It is important to note that the value of 
r]e is of order a few, i.e., the electron degeneracy is moderate. This justifies that the exact 
Fermi-Dirac distribution must be used in the calculations, and analytic approximations for 
either extremely degenerate electrons (r/e 3> 1) or fully nondegenerate electrons (r/e -C 1) are 
invalid for hyperaccretion disks. 

Figure 2 shows the composition of disk matter. In the outer region between R = 500Rg 
and R ~ 200Rg, almost all a-particles are not disintegrated, so that the a-particles fraction 
Ya keeps to be ~ 0.25, the electron fraction Y^ keeps to be ~ 0.5, and the free proton fraction 
Yp, the free neutron fraction Y^, and the free nucleon fraction Xnuc are all keeping ~ 0. From 
it! ~ 200Rg inwards, the disintegration of a-particles causes Y^, to decreases and Xnuc to 
increase dramatically. Because of the neutronization processes favored by high temperature, 
Ya greatly exceeds Yp, and Ye decreases accordingly. In the innermost region R < lORg, 
a-particles are almost fully disintegrated, i.e., ~ 0, X^-^c ~ 1; and the neutronization 
makes Y^ larger than 0.9, and Yp and Y^ smaller than 0.1. 
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Contributions to the total pressure P from free nucleons and ct-particles Pgas, from 
degenerate electrons and positrons Pg, from photon radiation Prad, and from neutrinos P^, 
are drawn in Figure 3. It is seen that Pe > -Pgas in the outer region R > lOOP^, and Pgas 
becomes dominant at smaller radii because both of the disintegration of a-particles and the 
decreasing degeneracy of electrons; Prad is of some importance only in the outer region, and 
Pi, is small even in the innermost region R < lORg where the disk becomes optically thick 
to neutrinos. 

The neutrino optical depth T,y. is plotted in Figure 4. Two comparisons are made in 
Figure 4a. First, the optical depth for electron neutrinos r,y_, is several times larger than 
that for yU-neutrinos r^,^ and r-neutrinos t,^^ , and only t,^^ can become larger than 2/3 in 
the innermost region R < lORg, therefore r^.^ ought to be taken as the representative of 
Tui- Second, the contribution to r^,^ from absorption Ta,!,^ is more important than that from 
scattering Ts^^,^ , this is a surprising result as the general understanding in the literature was 
that scattering off nucleons is a far more important opacity source than absorption(e.g., Di 
Matteo et al. 2002; Lee et al. 2005). We notice that in calculating the neutrino cooling 
rate due to the URCA processes guRCA (which is dominant over other neutrino cooling 
rates, see Fig. 4b) those authors used an approximate formula that is valid only in the 
nondegeneracy limit. In fact, the electron degeneracy causes guRCA to increase greatly as in 
our calculations. According to equation (46) of Kohri & Mineshige (2002), in the complete 
electron degeneracy limit Qvrca becomes extremely large because ^urca oc r)e^\ It is indeed 
the case as seen from Figure 4b that the contribution to Ta,ue from absorption relating to 
the URCA processes turca is most important, the contribution from neutrino annihilation 
into e~e~^ pairs Te-e+ is of some relative importance only in the very outer region, and the 
contributions from nucleon-nucleon bremsstrahlung Tbrcm and from the inverse process of 
plasmon decay Tpiasmon are totally negligible because the electron degeneracy is not very 
high. As to Ts^i^^ (Fig. 4c), the contribution through scattering off free neutrons r^^c is 
dominant over other contributions due to free protons Tp^,^^ , electrons Te^ve , and a-particles 
Ta,i^e {Ta,ue is too Small to be seen in the figure), obviously because of the richness of neutrons. 

Having r^,^ in mind, we now check the validity of our bridging formula for 1^, i.e., 
equation (42). In Figure 5, calculated using equation (42) is drawn by the solid line. Here 
M— 5Mq s~^is taken because only at such a high accretion rate the inner region of the disk 
can become noticeably optically thick. For comparison, the dashed line in the figure draws 
Ye in the case that the disk material is assumed to be globally opaque to neutrinos, which is 
calculated with (see Eq. [39]) 



and the dotted line draws Y^ if the material is globally transparent to neutrinos, obtained 




-1 



(46) 
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with the approximate formula of Beloborodov (2003), 



= ^(1 - ^nuc) + ^nuc[^ + 0.487( 



Q/2 - /Xe 



)]■ 



(47) 



These two equations do not include the effect of varying neutrino optical depth. It is clear 
that our equation (42) does represent a bridge between the optically very thick and optically 
very thin limits, and it provides a reasonable estimate of in the intermediate regime where 
Ye is overestimated by equation (46) and underestimated by equation (47). 

In Figure 6 we show the ratios of various coohng rates to the viscous heating rate Qvis- 
The photon radiation cooling is never important in hyperaccretion disks (it is practically zero 
when we try to calculate it). Advective cooling Qadv is dominant only in the outer region 
of the disk it! > 200it!g, because in this region the photodisintegration of a-particles has not 
started and the neutrino emission is weak. In the middle region 200i?g ^ -R ^ 50i?g, the 
a-particle photodisintegration cooling Qphotodis dominates, corresponding to a sharp decrease 
of Ya and increase of Xnuc(see Fig. 2). In the inner region R < 50Rg, the neutrino cooling 
Qi, becomes dominant as expected. 

We do not present here numerical results for the disk structure with varying values of 
a and M. Differences caused by changing these two parameters are only quantitative and 
have been discussed in, e.g., Kohri et al. (2005) and Chen & Beloborodov (2006). Briefly 
speaking, a smaller than 0.1 will make the density higher, the electron degeneracy higher, 
the electron fraction lower, and the neutrino-dominated region larger; and at still larger 
accretion rates, advection can become dominant over neutrino cooling again in a very small 
region {R < 5Rg for M—5Mq s~^, see Fig. 3 of Gu et al. 2006), because in that region the 
optical depth is very large and neutrinos are trapped in the disk. 



Having the neutrino cooling rate Qi,, the neutrino radiation luminosity (before annihi- 
lation) Li, is obtained as 



In our calculations the inner and outer edge of the disk are taken to be it!j„ = 3Rg and 
Rout — 500Rg, respectively. 

For the calculation of the neutrino annihilation luminosity we follow the approach in 
Ruffert et al. (1997), Popham et al. (1999), and Rosswog et al. (2003). The disk is modeled 



4. 



Neutrino Radiation and Annihilation Luminosities 




(48) 
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as a grid of cells in the equatorial plane. A cell k has its mean neutrino energy neutrino 
radiation luminosity l^. , and distance to a space point above (or below) the disk dk- The 
angle at which neutrinos from cell k encounter antineutrinos from another cell k' at that 
point is denoted as O^k'- Then the neutrino annihilation luminosity at that point is given by 
the summation over all pairs of cells, 

= E E I E + - COS e,,,f 

i k ^ k' ^' 

+ E^^.^ E I E - ^°«^^^'), (49) 

where = (l/127r2)[ao/c(mec2)'][(Cv.,,, - CA,uf + (Cy,., + Ca,.J'] , A^,^ - (l/67r2)(ao/c) 
{2Cy^^ — CljyJ, with Cy^y^ and CA,Ui given in §2.3.1. The total neutrino annihilation lumi- 
nosity is obtained by the integration over the whole space outside the black hole and the 
disk, 

/•oo poo 

^z.F = 47r / / l^TpRdRdZ. (50) 

JRg Jh 

Figure 7 shows Lj, (the thick dashed line) and L^-p (the thick solid line) with varying 
M (M = 3Mq, a = 0.1, and j = l.ScRg are kept). For comparison, these two luminosities 
calculated in our previous work (Gu ct al. 2006), where the electron degeneracy was not 
correctly considered and was taken to be equal to 0.5, are also given in the figure by the 
thin dashed line and thin solid line, respectively. It is clear that the electron degeneracy and 
the lower 1^ resulted from the neutronization processes indeed suppress the neutrino emission 
considerably, the resulting L^, and Li^-p arc reduced by a factor ~ 30% — 70% comparing with 
their overestimated values in Gu ct al. (2006). Even so, the correct L^^-p is still well above 
]^g50 gj,gg g-i p];ovidcd M ~ IMq s^"*^, and reaches to ~ 10^^ ergs s~^ when M ~ IOMq s~^. 
Therefore, based on the energy consideration, neutrino-cooled accretion disks can work as 
the central engine of GRBs. Note that our calculations are for a nonrotating black hole, 
both Popham et al. (1999) and Chen & Beloborodov (2006) have shown that a spinning 
(Kerr) black hole will enhance the neutrino radiation efficiency, this only strengthens our 
conclusion here. 

In Figure 8 the corresponding neutrino radiation efficiency r]iy{= Ly/ Mc^) and neutrino 
annihilation efficiency r]uv{= L^y/Ly ) are shown by the dashed line and solid line, respec- 
tively. It is seen that r]y does not change much with varying M; while r]yi> increases rapidly 
with increasing M, this is because at higher accretion rates more neutrinos are emitted and 
in turn have higher probabilities to encounter with each other. 
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To see the spatial distribution of neutrino annihilation luminosity, we plot in Figure 9 
contours of {2T:Rli,y) in units of (erg cm~^), i.e., the neutrino annihilation luminosity of a 
circle with cylindrical coordinates R and Z, for M= IMq s^^. The figure has two important 
implications. First, it demonstrates the strong focusing of neutrino annihilation towards the 
central region of space. In fact, by performing the integration of equation (50) out to each 
radius we get that nearly 60% of the total annihilation luminosity is ejected from the region 
R < 20Rg. Second, it shows that the annihilation luminosity varies more rapidly along the 
Z coordinate than along the R coordinate, indicating that this luminosity is anisotropic and 
most of the annihilation energy escapes outward along the angular momentum axis of the 
disk. 



5. Summciry and Discussion 

In this paper we wish to discuss whether the annihilation of neutrinos emitted from 
hyperaccretion disks can provide sufficient energy for GRBs, i.e., to estimate the neutrino 
annihilation luminosity Lj^p. To do this, we need to know the neutrino optical depth r^,., 
because it determines the neutrino cooling rate Qi, (eq. [9]) and the neutrino radiation 
luminosity L^, (eq. [48]). 

To calculate contributions to r^,. from various absorption and scattering processes of neu- 
trinos (eq. [16]), we need to know the composition and physical state of disk matter, namely 
the fractions of electrons Yc and free nucleons -'^nuc; and the electron chemical potential He 
(or electron degeneracy rjc), with given density p and temperature T. We give exact defini- 
tions of the free proton fraction Yp and free neutron fraction Y^^ and their relations to 1^ and 
Xjiuc; and get three equations (eqs. [42], [44], and [45]) from the conditions of /3-equilibrium, 
nuclear statistical equilibrium, and charge neutrality that describe the interdependence of 
Yp, Xyxuc^ Pci and t^^ . 

We prove that the electron degeneracy has important effects indeed, but mainly not 
in the sense that it suppresses the creation of neutrinos from e^e+ pairs and enlarges the 
electron pressure (Kohri & Mineshige 2002; Kohri et al. 2005). As seen from Figures 4b and 
3, the neutrino cooling due to e~e"^ pair annihilation qe-e+ is much smaller than that due to 
the URCA processes ^urca, and the electron pressure is important only in the outer region 
of the disk where the ct-particle disintegration has not started. Instead, the main effects of 
electron degeneracy are: (1) It increases Q'urca greatly, so that the corresponding absorption 
makes a major contribution to the neutrino optical depth (Fig. 4); (2) It along with the 
neutronization processes cause Y^ to become smaller than 0.1 in the inner region of the disk 
(Fig. 2) where the neutrino cooling is dominant (Fig. 6). 
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The resulting L^-p is considerably reduced comparing with that in previous works where 
the electron degeneracy was not considered and was taken to be 0.5, however it is still 
likely to be adequate for GRBs, and its spatial distribution is likely to be anisotropic. 

Perhaps the main limitation of our calculations here is that they are one-dimensional. 
As in most of previous works in the field, we do not study the vertical structure of the disk 
and instead use the simple 'one-zone' approximation of the vertically-averaged model. In 
particular, we do not consider the distribution and transport of neutrinos in the vertical 
direction of the disk, a problem that remains unsolved or even rarely touched (Sawyer 2003). 
Strictly speaking, a reliable quantitative evaluation of the neutrino annihilation luminosity 
should require two-dimensional calculations, in which the vertical structure and neutrino 
transport are self-consistently included. But in view of the fact that most popular models 
of normal accretion disks (the Shakura-Sunyaev disk model, the slim disk model, and the 
advection-dominated accretion flow model) are also one-dimensional and have been proved 
successful, our results here may provide a plausible, though rough, estimate of the luminosity 
of neutrino-cooled accretion disks. 
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Fig. 1. — (a) Density p, (b) Temperature T, and (c) electron degeneracy rjg as functions of 
radius R, with the black hole mass M = 3Mq, mass accretion rate M = IMq s~^, viscosity 
parameter a = 0.1, and the accreted specific angular momentum j = l.ScRg. 
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Fig. 2. — Electron fraction Y^, free proton fraction Yp, free neutron fraction Y^, a-particle 
fraction Y^, and free nucleon fraction Xnuc as functions of R, with the same constant param- 
eters as in Fig. 1. 
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Fig. 3. — Contributions to the total pressure P from free nucleons and a-particles Pgas, from 
degenerate electrons P^, from photon radiation Prad, ^-nd from neutrinos Pi, as functions of 
R, with the same constant parameters as in Fig. 1. 
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Fig. 4. — (a) Total optical depth for electron neutrinos t^^^ , total optical depth for r- 
neutrinos and yU-neutrinos t^^ ^ , absorption optical depth for electron neutrinos Ta,uc i S'-nd 
scattering optical depth for electron neutrinos r^, ; (b) Quantity Ta and its contributions 
from the URCA processes txjrca, from neutrino annihilation re-c+ , from nucleon-nucleon 
bremsstrahlung Tbrem, and from the inverse process of plasmon decay Tpiasmon; and (c) Quan- 
tity Ts,ue and its contributions due to free neutrons Tn,iy^ , due to free protons Tp^i,^ , due 
to electrons t^^,^^ , and due to ct-particles Ta,,^^ as functions of R, with the same constant 
parameters as in Fig. 1. 
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Fig. 5. — Solid, dashed, and dotted lines draw calculated with equations (42), (46), and 
(47), respectively, with the same constant parameters as in Fig. 1 except forM = 5Mq s~^. 
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Fig. 6. — Ratios of the advective cooling rate Qadvj a-particle photodisintegration cooling 
rate Qphotodis , and neutrino cooling rate Qy to viscous heating rate Qvis as functions of i?, 
with the same constant parameters as in Fig. 1. 
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Fig. 7. — Neutrino radiation luminosity L,^ (the thick dashed hne) and neutrino annihilation 
luminosity L^y (the thick solid line) for varying M . The overestimated L,^ (the thin dashed 
line) and L^y (the thin solid line) taken from Gu et al. (2006) are also given. The constant 
parameters M, a, and j are the same as in Fig. 1. 
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Fig. 8. — Neutrino radiation efficiency rji, and neutrino annifiilation efficiency Tyj^u for varying 
M , witfi tfie same M, a, and j as in Fig. 1. 
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R/R 

g 

Fig. 9. — Contours of the neutrino annihilation luminosity of a circle with cylindrical 
coordinates it! and Z. The number attaching to each line is this luminosity in units of 
(erg s~^ cm~^). The shaded region shows the accretion disk. The constant parameters are 
the same as in Fig. 1. 



